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TECHNICAL  NOTE  NO.  15 66 

DAMPING  IN  PITCH  AND  ROLL  OF  TRIANGULAR 
WINGS  AT  SUPERSONIC  SPEEDS 
By  Clinton  E.  Brown  and  Mac  C.  Adams 

SUMMARY 


A  method  is  derived  for  calculating  the  damping  coefficients  in 
pitch  and  roll  for  a  series  of  triangular  wings  and  a  restricted  series 
of  sweptback  wings  at  supersonic  speeds.  The  elementary  "supersonic 
source"  solution  of  the  linearized  equation  of  motion  is  used  to  find 
the  potential  function  of  a  line  of  doublets,  and  the  flows  are  obtained 
by  surface  distributions  of  these  doublet  lines.  The  damping  derivatives 
for  triangular  wings  are  found  to  be  a  function  of  the  ratio  of  the 
tangent  of  the  apex  angle  to  the  tangent  of  the  Mach  angle.  As  this 
ratio  becomes  equal  to  and  greater  than  1.0  for  triangular  wings,  the 
damping  derivatives,  in  pitch  and  in  roll,  become  constant.  The  damping 
derivative  in  roll  becomes  equal  to  one-half  the  value  calculated  for 
an  infinite  rectangular  wing,  and  the  damping  derivative  in  pitch  for 
pitching  about  the  apex  becomes  equal  to  3.375  times  that  of  an 
infinite  rectangular  wing. 


INTRODUCTION 


In  reference  1,  a  straightforward  method  was  found  for  calculating 
the  lift  and  the  drag  due  to  lift  of  triangular  wings.  The  present 
paper  extends  the  method  to  the  calculation  of  rolling  and  pitching 
motions  of  the  wings.  The  damping  coefficients  in  roll  and  pitch  for 
the  limiting  case  of  very  slender  wings  have  been  calculated  (reference  2). 
The  present  theory  is  not  limited  by  the  size  of  the  apex  angle,  and 
triangular  wings  with  leading  edges  ahead  of  and  behind  the  Mach  cone 
originating  at  the  apex  of  the  wing  are  treated. 

In  the  present  theory,  based  on  the  linearized  equations  of  motion, 
the  wing  is  represented  by  a  doublet  distribution  which  can  be  shown  to 
be  equivalent  to  a  vortex  distribution.  An  integral  equation  is  found 
which  can  be  easily  solved  by  analogy  with  known  relations  for  two- 
dimensional  incompressible  flow.  The  pressure  distributions  presented 
may  be  used  to  calculate  the  damping  coefficients  of  a  limited  series  of 
wings  for  which  the  trailing  edges  are  cut  off  so  that  they  lie  ahead  of 
the  Mach  cone  springing  from  their  foremost  point. 
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SYMBOIS 


x,  J,  z 
xl*  7i>  Z1 


p  =  /m2 

0 

00 

*D 

K 

A 

C 


m 


€ 

f  (a) 

•-i 

_  2, 


coordinates  of  field  point  (see  fig.  l) 
coordinates  of  a  source  or  doublet 

disturbance— potential  function 

potential  of  supersonic  source 

potential  of  supersonic  source  distribution 

potential  of  supersonic  doublet  distribution 

potential  of  a  line  of  doublets 

source  or  doublet  strength 

tangent  of  half— apex  angle 

jiff  force\ 


lift  coefficient 


|pT2s 


pitching-moment  coefficient 


'Pitching  moment'' 
|pV2Sc- 


^  .  /Rolling  moment 

rolling-moment  coefficient  /  - - — r - 

l  |pV2Sb 

half  of  apex  angle  of  wing 
doublet— line-distribution  function 


root  chord 

mean  aerodynamic  chord 
point  about  which  the  wing  pitches 


P  P 

(Local  chord)  dy  =  % 
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s 

a 

p 

t 

K 

w 

n 

E'  (pc) 


wing  area 

angular  velocity  of  pitch 
angular  velocity  of  roll 
maximum  span  of  wing 
constant 


z-component  of  velocity 

small  quantity 

complete  elliptic  integral 
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F*  ((3C) 


Subscripts: 

<1 

P 

i 


ANALYSIS 


Solutions  must  be  found  that  satisfy  the  linearized  differential 
equation  of  a  nonviscous  compressible  fluid  written 

o  (i 

ar2  as2 

where  x,  y,  z  are  Cartesian  coordinates  (see  fig.  l),  and  <f>  is  the 
disturbance-potential  function  created  by  the  wing.  An  elementary 
solution  of  this  equation  known  as  the  potential  of  a  supersonic  source 
may  be  written 


°  |/(x  -  xx)2  -  P2(y  -  yq)2  -  P2(z  -  zq)2 

The  quantity  A  is  the  strength  coefficient  of  the  source.  New 
solutions  may  be  obtained  by  superposition  of  such  potentials  as  shown 
in  reference  3.  For  example,  a  distribution  of  sources  oyer  a  portion 
of  the  xy— plane  would  give  the  potential 


X  -  Xq)^~— (32 (y  -  yqj^  -  P2z2 


(3) 


where  the  limits  chosen  must  be  such  that  all  sources  will  be  located 
within  the  forward  Mach  cone  from  the  field  point  (x,y,z).  Another 
solution  may  now  be  obtained  by  differentiation  with  respect  to  any  of 
the  coordinate  directions,  that  is. 
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Sz 


cLXi 

|/(x  -  x-l)2  -  (3 2(y  -  y-^)2  -  p2z2 


(4) 


This  solution^  however ,  may  "be  considered  "the  vertical  or  z— component 
velocity  of  the  source-distribution  potential  and  as  shown  in 

reference  3 


±*A(x,y) 


(5) 


The  step  taken  in  equation  (4)  also  corresponds  to  the  formation  of  a 
doublet  potential,  that  is,  ^  represents  a  distribution  of  doublets 
over  the  xy-plane  with  strengths  proportional  to  A(x1,y1).  Tor  any 
known  doublet  distribution,  the  velocity  component  parallel  to  the  surface 
in  any  direction  s  may  immediately  be  obtained  from  equation  (5) 


±z 


(6) 


The  foregoing  results  are  analogous  to  incompressible— flow  relations  and 
it  may  be  stated  in  general  that  for  every  doublet  distribution  there  is 
a  vortex  distribution  which  will  produce  a  similar  flow.  The  vortex 
distribution  and  doublet  distribution  are  directly  related  by  equations  (5) 
and  (6).  These  simple  concepts,  given  first  by  Prandtl  (reference  4),  may 
be  used  directly  to  obtain  the  solution  of  problems  in  which  the  pressure 
distributions  are  given,  such  as  airfoils  of  uniform  loading.  If  the 
equation  of  the  surface  is  given  and  the  pressure  distribution  is  required, 
integral  equations  must  be  solved.  In  certain  cases,  the  problem  may  be 
simplified  if  the  form  of  the  final  potential  is  known.  In  reference  2 
the  disturbance  potential  for  wings  of  very  low  aspect  ratio  was  found  to 
be  in  the  form 


*  *  A(y) 


(T) 
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This  form  of  the  potential  appears  quite  logical  from  the  standpoint  of 
satisfying  the  boundary  conditions  for  steady  rolling  or  pitching.  In 
the  following  analysis,  the  assumption  of  a  potential  in  the  form  of 
equation  (7)  is  shown  to  be  correct;  however,  it  should  be  pointed  out 
that  the  potential  of  this  type  must  be  restricted  to  the  linearized 
theory  and  is  not  of  the  3ame  general  nature  as  that  of  a  conical  field 
which  exists  even  in  the  nonlinear  problems. 

From  equation  (7)  the  doublet  distribution  over  the  surface  will 
be  in  the  form 


(8) 


and  under  the  assumptions  of  the  linearized  theory  the  lifting-pressure 
coefficient  is  now: 


The  formation  of  the  integral  equation  follows  the  method  of 
reference  1.  A  potential  that  represents  a  line  of  doublets  in  the 
xy-plane  at  an  angle  tanker  to  the  x-axis  is  derived  in  the  form  of 
equation  (7)-  Use  is  made  of  the  boundary  conditions  to  set  up  an 
integral  equation  that  introduces  the  unknown  distribution  function  f(cr). 

The  potential  of  the  doublet  line  may  be  obtained  by  following  a  procedure 
similar  to  that  used  in  obtaining  equations  (3)  and  (4),  and  by  substituting 
the  expression  for  A  given  in  equation  (8)  into  equation  (4).  The 
expression  obtained  in  the  following  equation  may  be  seen  to  represent  a 
line  of  doublets  along  which  the  doublet  strength  increases  as  x3: 


-1  r  ~xi3 

^0  i/(x  -  xi  )2  -  p2(y  -  ot/l)2'  -  02z2 


2z(x  ~  P2cry )  f 

-P2o2)5/2\^ 


+  2p2z  y/x 2  —  p2(y2  +  z2) 

(1  -  pSo2)2 


(10) 
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where 


_ (x  ~  P2cry) _ 

[/l  —  P2^2  /x2  -  P2(y2  +  z2) 


and  x1  is  the  value  of  Xj_  for  which  the  denominator  of  the  integrand 
vanishes.  The  potential  of  the  complete  wing  may  now  he  obtained  by  an 
integration  with  respect  to  the  dimensionless  parameter  a 


0 


(n) 


where  tan- ^-C  =  €,  the  half -apex  angle,  and  f(cr)  is  an  unknown  distri¬ 
bution  function.  The  z— component  velocity  w  can  be  written  for  P  ^ 

approaching  zero 


Pf(g)(i  -  P2CT0) 

(1  -  32o2)5/2 


3  coth  H - 5-^ - 

?2  -1 


+ 


2x 


pc  pf  (cr)  1/1  -  p2e2 
PC  (1  -  p2cr2)2 


d(p<r) 


(12) 


where  0  =  ^  for  convenience.  The  boundary  conditions  for  rolling  may 
now  be  written 


w  =  -py 


or 


w 

X 


=  -p0 


(13) 


For  pitching  about  the  y-axis,  there  is  obtained 
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or 


(14) 


Introduction  of  equations  (13)  and  (l4)  into  equation  (12)  provides 
integral  equations  which  theoretically  can  "be  solved  for  the  unknown 
function  f(a).  Simpler  relations,  however,  may  he  obtained  if 
equation  (12)  is  differentiated  twice  with  respect  to  0  to  obtain  the 


quantity  £lsAl. 

bed 


The  method  for  differentiating  is  indicated  in  the 


appendix  and  gives 


d2(w/x) 


lim 


6 


/ 


2  2 
i  -  pe 


r ^e~r>)  p3f(g)  A(Pq) 

-PC  (per  -  pe)4 


6f L  -  p 


2q2 


^PPC  p3f (cr)  d(pa) 
p(e+T))  (pa-pe)4 


-  kj±~- 


p2©2 


f"(0)  f(0) 


(15) 


The  "boundary  conditions  require  the  foregoing  quantity  to  be  zero  for 
both  rolling  and  pitching  with  the  additional  requirements  on  f  (cr)  that, 
for  rolling,  at  the  point  0=0 


(w/x)p  =  0 


(16) 


and,  for  pitching. 


d(w/x) 

be 


=  o 


(17) 


Equation  (15 )  now  yields,  for  rolling. 


lim 

T)-^0 


(0-ti)  f (cr)-  do 


6 


C  f(oL  do 


(o  -  0) 


S  +  u  J,.  ^  : - T7£ 


( 0+T]  )  (o  -  0)' 


f”(e)p  f(0)p 


=  o  (18) 
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and,  for  pitching. 


lim 

TJ - >0 


f  r(6- t))  f(tf)a  dcr 

/  . <1  .  ,  £ 

pC  f(cr)  da 

1  I. 

f"(e)q  f(e)q]] 

'■M3  (a  -  0)^  ^ 

(0+n)  (cr-0)1^ 

i.  *.>Ji 

=  0  (19) 


Equations  (l8)  and  (19)  are  identical  to  the  equations  that  would  he 
obtained  for  similar  boundary  conditions  on  a  two-dimensional  flat  plate 
if  an  analogous  process  of  distributing  the  doublets  were  followed.  (See 
appendix.)  The  analogue  for  the  rolling  motion  of  a  triangular  wing 
would  be  a  two-dimensional  flat  plate  rotating  about  its  midchord  point 
in  a  stationary  stream.  The  surface  potential  distribution  and  therefore 
the  doublet  distribution  would  be 


f(c)p  =  Kpff  l/S2"-  ct2 


(20) 


For  the  pitching  condition  the  analogue  would  be  a  two-dimensional 
flat  plate  in  a  stream  flowing  normal  to  the  surface.  The  potential  or 
doublet  distribution  would  be 


*<*>q 


^ 


(21) 


These  potentials,  which  can  he  found  in  references  2  and  5,  satisfy 
equations  (l8)  and  (19)  hy  analogy;  however,  the  conditions  of  equations  (l6) 
and  (17)  must  he  shown  to  he  satisfied.  For  the  calculations  of  (w/xL 
8(w/x)  v 

and  — ^ - ,  and  the  evaluation  of  Kp  and  K^,  only  one  value  of  0 


need  he  considered.  This  value  may  conveniently  he  set  equal  to  zero.  Far 
rolling  motion,  equation  (20)  indicates  the  doublet  distribution  to  be 
antisymmetric.  Therefore  the  value  of  w/x  at  0  =  0  must  he  zero,  and 
the  condition  of  equation  (l6)  is  satisfied.  For  the  pitching  motion, 
the  doublet  distribution  is  symmetrical  about  6=0  and  therefore  the 


d0 

is  satisfied. 


quantity 


must  be  zero  at 


0=0  and  the  condition  of  equation  (17) 


The  constants  K  and 
P 

obtained  in  the  appendix  for 


may  now  be  evaluated  from  the  relations 
6  =  0 


mi*: 
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a(-w/x) 

de 


+  SK, 


f 


PC  /pgC2  -  P2^ 
-pc  (l  -  P2a2)2 


d(p<x)  +  K, 


I 


P(0-^1 ) 

-pc 


l/p2c2  -  pV 

— - —  d(Pcr) 

pVCl  -  P2a2)2 


f'PG  |/p2C2  -  P^cr2 
+  pV(l  - 


d(Pc)  - 


gpCKq 

Ptj 


(22) 


=  -P  = 


lim 


^  Cl  tai*"1 1/1 " p2a2  d(pa) 


r3(0^l)  i/q2q2  _  e2a2  N  rPC  l/p2C2  -  p2a2 

22^  /  -^-2 - 2_2_  d(Pcr)  +  2K  /  - — - - d(pa) 

^J-£C  p2cj2  JP(e+ii)  P2^2 


rpc  p2q2/p2c2  -  P2,g£  d(Pa)  _  k_  fpc  v^gc2  -  pjg.  d(3ff) 

-pc  _  p 2ar2^2  “PC  (1  -  P2^ 
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Equations  (22)  and  (23)  may  be  Integrated  by  use  of  tables  (reference  6) 
to  give 


P  =  «Kp 


2  -  p2C2 
1  -  p2C2 


E*  (pc) 


P2C2 

1  -  p2c2 


F*  (pc) 


(2k) 


q  =  jtE. 


2  2 

1  -  2p  C^ 


1  -  P2C2 


C*(pC)  + 


2  2 
P  cd 


1  -  p2c2 


F*  (pc) 


(25) 


F1 (PC)  and  E»(PC)  are  complete  elliptic  integrals  of  the  first  and 
second  kind. 


The  pressure  distribution  for  the  rolling  wing  may  now  he  obtained 
from  equations  (9),  (20 ),  and  (24)  and  the  pressure  coefficient  is 


P  = 


kxpc2e 


2  -  P2^ 


1  -  p2c2 


E*  (pc)  - 


P2C2 


i  -  p2e£ 


F‘(pC) 


Jc2  -  e2 


(26) 


Integration  of  the  pressures  over  the  wing  surface  gives  the  forces  and 
moments  acting  on  the  wing.  The  nondimens ional  derivative  may  then 

be  found  P. 


C2  = 


— rcC 


2  -  p2C2 


P2C2 


E*  (pc)  - 


P2C2 


(27) 


2  2 

1  - 


F*  (pc) 


In  the  analysis  the  pitching  axis  has  been  taken  at  the  wing  apex; 
however,  in  application  it  is  desirable  to  obtain  the  pressure  distribu¬ 
tion  and  the  force  and  moment  coefficients  for  pitching  about  any  point. 
A  superposition  of  motions  is  therefore  required.  The  pitching  motion 
about  any  point  xQ  can  be  made  up  of  a  pure  pitching  motion  about  the 
apex  of  the  wing  combined  with  a  vertical  translational  motion  of 
velocity  qxQ.  The  pressure  distribution  for  this  translational  motion 

corresponds  to  that  of  a  wing  at  a  constant  angle  of  attack  of 


22 
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q*o 


(See  references  1  and  7-) 

<1*0 


constant  angle  of  attack 


is 


The  pressure  distribution  for  the 


P  = 


-4C2gXQ 


ye* (pc)  i/cr  -  e2 


(28) 


Combining  equations  (9)>  (21 ) ,  (25),  and  (28)  gives  for  the  pressure 
distribution  in  the  pitching  case 


P  = 


4qx 


v  /c2  -  e2 


2C2  _  q2 


*0C‘ 


1  -  202C2 


2  2 
i  -  rc* 


E»(pC)  + 


P2C2 


P2c2 


F*  (pC) 


xE* (PC) 


(29) 


Integration  of  the  pressures  over  the  wing  surface  and  formation  of  the 
nondimensional  derivative  yields 


6*C 


l-2|gj£E»(pc)  + 

1  -  P^CT 


p2c2 


F*  (PC) 


^jcCxq 

Ef(pC)c 


and 


1  -  2P2C2 

1  -  p2c2 


E* (pc)  + 


1  -  p2c2 


F*  (pc) 


+ 


4nCxQ^l  - 
cE’ (pc) 


(30) 


(31) 


where  c"  is  the  mean  aerodynamic  chord. 

Calcinations  of  these  derivatives  for  triangular  wings  having  their 
leading  edges  outside  the  Mach  cone  are  most  easily  made  by  the  source 
distribution  method.  In  this  method,  the  upper  and  lower  sides  of  the 
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wing  may  be  considered  independent  of  each  other.  The  source  distribution 
function  for  the  rolling  wing  is 


Sp 


(32) 


whereas  that  for  the  pitching  wing  is 


Sg 


cc 


Kx- 


(33) 


The  calculation  of  the  pressure  distribution  is  not  presented,  since  the 
subject  of  the  integration  of  source  distributions  has  been  well  covered 
in  reference  3* 

The  pressure  distribution  for  rolling  wings  outside  the  Mach  cone  has 
been  calculated  to  be 


P  = 


4pC2x 


*y(p2c2  -  1) 


372 


(1  +  p2C0) 


-1  1  +  p2c@ 

COS  - 

P(C  +  0) 


-  (1  -  p2ce) 


_i  1  -  p2ce 

COS  - 

p(c  -  0) 


(34) 


Integrating  the  pressures  over  the  wing  and  expressing  the  derivative  in 
nondimen3ional  form  gives 


C 


(35) 


For  the  pressure  distribution  due  to  pitching  about  the  point  Xq,  a 
combination  of  flow  patterns  must  again  be  used.  The  pressure  distribution 

qxQ 

of  a  wing  at  uniform  angle  of  attack  is  (reference  3) 


P  = 


4<i^C 


ff/\ 32C2 


-  1  L 


COS 


1  1 


p2C0 


-1  1  +  P  C0 
+  cos  - e - 

P(c  -  0)  p(c  +  0) 


(36) 


The  pressure  distribution  for  pitching  then  becomes 
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P  = 


4qx 

rf\Tp 


2pci/4  -  p2e2 


p2c2 


p3c3  -  2pc  -  p0 

( p2c2  -  1)3/2 


COB 


1 


1  +  pgce 
P(c  +  0) 


p3c3  _  2 pc  +  pe  -1  1  -  p2ce 

+  - — ; COS  — : r  1  -r 

(P2C2-  l)3/2  3(C-0) 


4qx0pc 

TtVp/p2C2  -  1 


1  +  p2ce 
P(c  +  0) 


+  cos 


-1 


1  -  p2ce 

P(c  -  0) 


(37) 


The  nondimens ional  derivatives 


CT  and  Cm  then  become 

q 


8  8xo 

P  pc 


(38) 


P  Pc 


(39) 


DISCUSSION  AND  CONCLUSIONS 


Expressions  for  the  lifting-pressure  coefficients  over  triangular 
wings  in  roll  are  given  in  equations  (26)  and  (3*0  and.  in  pitch  in 
equations  (29)  and  (37)-  Equations  (2 6)  and  (29)  are  for  wings, 
inside  the  Mach  cone  and  equations  (3*0  and  (37)  for  wings  outside 
the  Mach  cone-  Typical  pressure  distributions  are  shown  in  figure  2 
in  which  the  pressure  distributions  for  the  two  wings  in  pitch  are 
for  pitching  about  the  apex- 

Expressions  for  the  quantities  C^,  Cj^,  and  C^^  are  given  in 

equations  (27),  (30),  and  (3l),  respectively,  for  the  case  of  the  wing 
inside  the  Mach  cone  and  in  equations  (35),  (38),  and  (39)  for  wings  lying 
outside  the  Mach  cone.  It  will  be  seen  that  the  parameters  PC^,  PCL^, 

and  PC™  may  be  expressed  as  functions  of  PC  where 

m.-, 
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PC 


tan  e 
tan  n 


The  stability  derivatives  may  therefore  be  plotted  against  this  parameter, 
to  give  curves  which  will  hold  for  all  triangular  wings  at  any  Mach  number. 
These  curves  are  given  in  figures  3,  and.  5.  For  values  of  PC 
approaching  zero  the  values  of  the  derivatives  closely  approach  those 
given  in  reference  2  which  were  based  on  the  assumption  of  very  low 
aspect  ratio. 


For  values  of  pc  ^  1  (that  is,  for  the  wing  lying  outside  the  Mach 

cone),  the  quantities  0C 1  and  PC^  become  constant,  and  equal  to 

P  q 

—  and  — 1,  respectively,  (the  pitching  being  about  the  2-c  point).  In 
comparison,  the  value  of  PC}  and  PCL  for  infinite-span,  rectangular 

o  8  ?  ^ 

wings  are  -  j  and  —  — ,  respectively,  (the  pitching  being  about  the 
leading  edges). 


It  should  be  pointed  out  that  the  pressure  distributions  given  in 
this  paper  may  be  used  directly  to  calculate  the  damping  in  pitch  and 
roll  for  wings  having  trailing  edges  cut  off  ahead  of  the  Mach  cone,  the 
most  interesting  of  this  series  being  the  so-called  "arrow  wings." 


It  is  apparent  that  a  suction  force  exists  at  the 
wings  in  pitch  and  roll  whenever  the  leading  edges  are 
Mach  cone.  A  method  for  obtaining  the  values  of  these 
was  derived  in  reference  1. 


leading  edges  of 
swept  behind  the 
suction  forces 
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Introducing  the  limits  and  then  setting  z  =  0  gives 
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Equation  (12)  may  now  be  rewritten  for  w/x  with 


Differentiating  twice  with  respect  to  j  gives 
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Figure  1.-  Coordinate  system. 


Figure  2.-  Pressure  distributions  for  rolling  and  pitching  about  apex. 
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Figure  4.-  Stability  derivative  Cr 


about  the  ^-c  point  for  triangular  wings. 


